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From an analysis of various types of data obtained in rel- 
ativistic nuclear collisions, the following picture has emerged 
in thermal and hydrodynamical descriptions: as the fluid ex- 
pands and cools, particles first undergo a chemical freeze- 
out at T c h.f. ~ 160 — 200 MeV then a thermal freeze-out 
at T t h.f. ~ 100 — 140 MeV. In this paper we show how to 
incorporate these separate freeze-outs consistently in a hy- 
drodynamical code via a modified equation of state (general 
case) or via a modified Cooper-Frye formula (particular case 
of T c h.f. close to Tth.f. or few particle species undergoing early 
chemical freeze-out). The modified equation of state causes 
faster cooling and may have sizable impact on the predicted 
values of observables. 
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I. INTRODUCTION 

The behavior of strongly interacting matter under ex- 
treme conditions of pressure and temperature is the 
subject of the research programs at CERN (SPS) and 
Brookhaven (AGS and RHIC). Hydrodynamical and 
thermal models have been used extensively to describe 
data from these collisions and the following picture has 
emerged (see e.g. from a study of collisions at SIS, 

AGS, SPS energies, with a variety of targets and projec- 
tiles (for a compilation see e.g. ||). It is also expected 
to hold at RHIC and LHC energies. The dense and hot 
fluid expands and cools until chemical freeze-out occurs 
for some species of particles. Namely these particles stop 
having inelastic collisions so that their abundances are 
frozen. Therefore by studying the abundances of these 
chemically frozen particle species, one learn the condi- 
tions at chemical freeze-out. For example, at CERN, 
T ch . f . ~ 160 - 200 MeV [@. The fluid goes on cool- 
ing until thermal freeze-out happens. Precisely, particles 
stop having elastic interactions and so the shape of their 
momentum distribution is fixed. Therefore by studying 
these spectra, one gets information about the conditions 
at thermal freeze-out. For example at CERN, using vari- 
ous types of particles || or combining information about 
the spectrum of a single particle species and its Bosc- 
Einstein correlations one extracts Tth.f. ~ 100 — 140 
MeV. There are some deviations to this picture. For ex- 



ample, some particles like the Q may undergo both freeze- 
outs almost together and early (due to their small cross 
section). 

Though this picture is simple and consistent with data, 
its theoretical justification needs further scrutiny. First, 
in a way, this picture works too well as a statistical de- 
scription seems to apply even to very elementary systems 
||. This has been debated a lot yet remains an open 
question. Second, in this description, it is assumed that 
particles make sudden freeze-outs. For example, when 
they cross the 180 MeV temperature three-dimensional 
surface in the fluid they immediately stop interacting in- 
elastically. In reality one expects that this should happen 
over a certain length. A formalism to account for finite 
freeze-out volumes and the subsequent new data interpre- 
tation have been presented in || . Third, in the freeze-out 
scenario, when computing particle distribution with the 
usual Cooper-Frye formula fTc[, there may be particles 
contributing negatively, corresponding to particles that 
are in the frozen out region and re-entering the interact- 
ing region. Physically this should not happen but in the 
calculation, it may. Ways to deal with this problem and 
how it affects particle data interpretation can be found 
ml!. 

Leaving aside these problems for further studies, the 
objective of this paper is the following. In the past, 
cascade event generators (Fritiof, Venus, RQMD, ARC, 
etc) were employed in particular by experimental groups, 
to study their data. Recently, simple thermal and 
hydrodynamics-inspired models have been used increas- 
ingly. It seems useful therefore to start developing more 
sophisticated hydrodynamical codes to extract phys- 
ical information from data. Though data call for sepa- 
rate chemical and thermal freeze outs, no hydrodynami- 
cal code so far includes chemical and thermal freeze outs 
self-consistently. Namely, the effect of the early chemi- 
cal freeze out on the fluid expansion is never taken into 
account. In this paper, we discuss how to incorporate 
separate chemical and thermal freeze-outs in a hydrody- 
namical code and show that in certain cases, this will 
have sizable impact on the predicted values of observ- 
ables. 

Finally let us mention that early universe and rela- 
tivistic heavy ion collisions call for different treatments 
for the following reason. In the early universe, the ex- 
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pansion rate (H ~ 10 4 s~ 1 ) at the QCD phase transi- 
tion is many orders of magnitude smaller than for rel- 
ativistic heavy ion collisions (H ~ 10 21 — 10 23 s _1 ). 
In both typical reaction rate for hadrons is 

r = <rnv re i ~ XfmT 1 — 10 s . Therefore while in 
the early universe, chemical and thermal equilibrium be- 
tweeen hadrons must have prevailed (r >> H), this is 
not the case for relativistic heavy ion collisions (T may 
be ~ H). Moreover in relativistic heavy ion collisions, 
we expect separate chemical and thermal freeze-outs (i.e. 
at distinct temperatures) because inelastic collisions (re- 
sponsible for chemical equilibrium) need higher center of 
mass energy to be operative in general, compared to elas- 
tic collisions (responsible for thermal equilibrium). This 
distinction needs not be made in the early universe (both 
rates Ti ne i and r e ; being much larger than the expansion 
rate). 



II. INCLUSION OF SEPARATE CHEMICAL AND 
THERMAL FREEZE-OUTS IN THE 
HYDRODYNAMICAL EQUATIONS 

To simplify the discussion, we use a known and simple 
hydrodynamical model, Bjorken one-dimensional boost 
invariant model jD| . Before chemical freeze-out, the 
fluid evolution is governed by the hydrodynamical equa- 
tions 

dt e + p 
dt + ~T 
dn B , n B 



= 



at ' t = ' 

The last equation can be solved easily 

nB{to)to 
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The first equation must be completed by the choice of an 
equation of state p(n B , e), for the pressure as function of 
the net baryon density and energy density. 

Also to simplify the discussion, we suppose that chem- 
ical or thermal freeze-out occurs at some fixed temper- 
ature (as often assumed in the analysis of experimental 
data). Attempts to incorporate more physical freeze-out 
conditions have been carried out fl^- [l7| , [l2]| and in princi- 
ple might be incorporated in the scheme described below. 

When the fluid temperature has decreased to some 
temperature T c h.f., (which corresponds to a certain time 
tch.f.), some particle species get their abundances frozen. 
To fix ideas, we suppose that A and A are in this sit- 
uation. Then in addition to the above hydrodynamical 
equations, we introduce separate conservation laws for 
these two types of particles for time t > £ c w., namely 

dn A t n A n 



dt 
dt 



t 

2* = o. 



(5) 



These equations have solutions of the same form as (3) 
but with to substituted by t c h- Therefore what remains 
to be done is to solve the energy-momentum equation 
(1) with a modified equation of state, to account for the 
particles who make an early chemical freeze-out. 

We suppose that the fluid is a gas of non-interacting 
resonances. Then for particle species i, using an expan- 
sion in term of modified Bessel functions pl| to allow the 
study of their limit more easily in (9-11), 
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where m% is the particle mass, gi, its degeneracy and /Uj, 
its chemical potential, the minus sign holds for fcrmions 
and plus for bosons. In principle each particle species i 
making early chemical freeze-out has a chemical poten- 
tial associated to it; this potential controls the conser- 
vation of the number of particles of type i. For particle 
species not making early chemical freeze-out, the chem- 
ical potential is of the usual type, yn = Bifi B + SiHs, 
where fi B (ps) ensures the conservation of baryon num- 
ber (strangeness) and Bi (Si) is the baryon (strangeness) 
number of particle of type i. So the modified equation of 
state depends not only on T and \x B but also /ia, etc. 
(the notation "etc" stands for all the other particles mak- 
ing early chemical freeze-out) [[l9| . This complicates the 
hydrodynamical problem, however we can note the fol- 
lowing. 

If m,i — fii » T (the density of type i particle is low) 
and rrii >> T, (these relations should hold for all parti- 
cles except pions and we checked them for various times 
and particle types) the following approximations can be 
used 
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Pi are written in term of n, 
can work with the variables 



We note that and 
and T. Therefore we 
T, n B ,n A ,n A ,etc, rather than T, [i B , fj,\, fi A , etc. The 
time dependence of n\,n A , etc is known as discussed al- 
ready. So the modified equation of state can be computed 
from t, T and fi B . 

The scheme presented above can easily be generalized 
to particles making chemical freeze-out at different times 
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(using different t c h.f.) and particles doing chemical and 
thermal freeze-outs together (whose contribution drop 
out of the equation of state). One can show that entropy 
is conserved even in the presence of an early chemical 
freeze-out when the hydrodynamical equations are satis- 
fied by a perfect fluid [|l9| . 

For illustration, we present results using in the equa- 
tion of state, the basic multiplets of resonances (pseu- 
doscalar meson octet plus singlet, vector meson octet 
plus singlet, baryon octet and baryon decuplet) and sup- 
posing that the early chemical freeze-out occurs at 180 
MeV, a value typical for SPS energy more or less indepen- 
dently of the projectile ||. We use for initial conditions 
T = Mso = 200 MeV, r = 1 fm, so that fi B ch.f. = 210 
MeV, which is in agreement with results for S or Pb at 
SPS ||. In figure 1, we compare the behavior of T and 
/iB as function of t, obtained from the hydrodynamical 
equations using the modified equation of state and the 
unmodified one. For the modified equation of state, we 
considered two scenarios: (I) all strange particle in the 
basic multiplets, (II) all strange particles except K and 
if *'s, make an early chemical freeze-out. This is a conser- 
vative estimate, it is possible for example that the pions 
make an early freeze-out [ p"2"| . Comparing scenarios I and 
II, we note that if more particles undergo early chemical 
freeze-out, stronger effects for T(t) and /ze(t) are seen. 
We concentrate on I thereon. We see that the deviations 
between I and the unmodified equation of state case in- 
crease with time. In particular from this figure, if the 
thermal freeze-out occurs at 110 MeV, the thermal freeze- 
out time is 13 fm for the modified equation of state and 
20 fm for the unmodified one; the corresponding bary- 
onic potentials are not very different, 405 and 375 MeV 
respectively. If the thermal freeze-out occurs earlier, say 
at 140 MeV, the difference in the thermal freeze-out times 
would be much less. An immediate consequence of this, 
is that the thermal freeze-out volume (in our case sim- 
ply proportional to time) may be much smaller for the 
modified equation of state. We expect this conclusion to 
hold qualitatively even in the presence of transverse ex- 
pansion: in this case, expansion is faster and the thermal 
freeze-out temperature is reached faster, so the thermal 
freeze-out times for the modified and unmodified equa- 
tion of state are less different, however there is a compet- 
ing effect for the thermal freeze-out volumes, they now 
scale with higher powers of time. We also expect this 
conclusion to hold at RHIC energies, with T c h still of or- 
der 180 MeV (it cannot be much higher since a transition 
to quark-gluon plasma is expected at about this temper- 
ature from lattice gauge simulations) but a lower value 
°f I^B.ch (less baryon stopping is expected at RHIC than 
SPS), the value of T t h may be a little smaller than at SPS 
pol , [l5| and ^s,th will be higher than pLB,ch (cf. figure 1). 
We conclude that if the chemical and thermal freeze-out 
temperatures are very different (in our simplified case, 
180 and 110 MeV) or if many particle species make an 



early chemical freeze out (e.g. also pions), it is impor- 
tant to take into account the effect of the early chemical 
freeze- out on the equation of state to make predictions for 
observables which depend on thermal freeze-out volumes, 
for example particle abundances and eventually particle 
correlations. This is our main result. 
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FIG. 1. hb and T as function of time in the case where all parti- 
cles have simultaneous freeze-outs (dashed line) and (I) all strange 
particles in basic multiplets make an early chemical freeze-out (con- 
tinuous line), (II) all strange particles except K and K*'s make an 
early chemical freeze-out (dotted line). 



III. PARTICULAR CASE WHERE T CH ~ T TH OR 
FEW PARTICLE SPECIES UNDERGO EARLY 
CHEMICAL FREEZE OUT 

If the chemical and thermal freeze-out temperatures 
are not very different (say 180 and 140 MeV), or if few 
particle species make the early freeze out, one can pro- 
ceed as follows. One can use an unmodified equation of 
state in a hydrodynamical code and to account for early 
chemical freeze-out of species i, when the number of type 
i particles was fixed, use a modified Cooper- Frye formula 



En _ Ni(T ch . f .) 



dp 3 



Ni(T th .f.) Js th . f 



da^f(x,p). (12) 



The second factor on the right hand side is the usual 
one and it gives the shape of the spectrum at thermal 
freeze-out, the first factor is a normalizing term intro- 
duced such that upon integration on momentum p, the 
number of particles of type i is Ni(T c h.f)- For illustra- 
tion, we show results obtained with the hydrodynamical 
model, HYLANDER-PLUS @. It provides a numerical 
solution of the relativistic hydrodynamical equations in 
(3+1) dimensions with axial symmetry (for details, see 
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pl| p4[). It gives a good description of single-particle- 
rapidity data, transverse- momentum spectra of h~ , 7r~, 
p, p, K°, 7r~ /tt + and pion correlation data at CERN en- 
ergies, for an appropriate choice of the initial conditions, 
an equation of state incorporating a first order phase 
transition and a freeze-out temperature of 139 MeV. In 
figures 2 and 3, we show results obtained for A, A and 
S, S respectively (neglecting resonance decays) and data 
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FIG. 2. Results from HYLANDER-PLUS for the transverse 
momentum spectra of A and A compared with data [22] for si- 
multaneous freeze-outs at T c h.f. = Tth.f. = 139 MeV (dash-dotted 
line) as well as separate freeze-outs at T c ^ f = 176 MeV (contin- 
uous line) or T ch f ^ = 184 MeV (dashed line) and T th f ^ = 139 
MeV. 
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FIG. 3. Same as figure 2 but for S and S. 

We see that both their shapes and abundances can be 
reproduced for T chJ , = 176 MeV and T fhJ , = 139 MeV, 
while simultaneous freeze-outs at T c/l j. = T t h.f. = 139 
MeV would yield the correct shapes but too few parti- 
cles. Therefore,results with HYLANDER-PLUS and a 
modified Cooper-Frye formula (12) support the separate 
freeze-outs picture. However in this code, Tth.f. is fixed 
to 139 MeV while as already mentioned, some data seem 
to imply lower thermal freeze-out temperatures. In the 
next generation of hydrodynamical codes, it is desirable 
to consider a wider range of T t h.f.- 

IV. CONCLUSION 



In summary, we showed how to incorporate separate 
chemical and thermal freeze-outs in a hydrodynamical 
code via a modified equation of state in section II (general 
case) or via a modified Cooper-Frye formula in section 
III (particular case oiT c h.f. close to Tth.f. or few particle 
species undergoing early chemical freeze-out). The mod- 
ified equation of state causes faster cooling and may have 
sizable impact on the predicted values of observables. 
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